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inferred, where the main contribution to the spectral line shape originates from the interaction with 
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I. INTRODUCTION 

Spectroscopy is one of the most powerful methods in plasma diagnostics. In particular, the analysis of the shape 
of spectral lines allows to determine the properties of the plasma such as temperature, density, and composition. For 
dense, strongly coupled plasmas, optical spectra have been investigated to infer the parameter values not only of 
laboratory, but also of astrophysical plasmas, see, e.g. Refs. [H, H, and Refs. respectively. 

The shape of spectral lines is determined by different processes. Besides the natural line width given by the finite 
lifetime of excited states, Doppler broadening is related to the thermal motion of the emitters. The influence of the 
surrounding medium becomes more important at increasing density and leads to pressure broadening. A large effect is 
caused by charged particles in the plasma. This Stark broadening is determined by the distribution of ions as sources 
of the microfield, whereas the contribution of the electrons is usually treated in impact approximation [1, Q . The 
Weisskopf radius can be used to distinguish between weak and strong collisions. However, also the neutral components 
in the surrounding medium contribute to line profile, i.e. as van der Waals broadening. 

A systematic approach to the shape of spectral lines in dense plasmas has been worked out on the basis of pertur- 
bation theory [l(|, and many-particle effects have been incorporated. Griem [ll[ gives a review of different classical 
and semiclassical methods as well as their application to line shapes in different plasmas. The unified theory [l2l . [l3j 
was developed to describe the center as well as the wings of the lines. A quantum statistical approach based on the 
many-particle Green functions method was presented by Giinter 0, [l5[ and has been developed systematically for 
charged perturbers. Special attention was payed to show the equivalence between this rigorous approach and the 
aforementioned classical and semiclassical concepts in the appro pria te limits. The quantum statistical approach has 
been applied to a broad range of systems, e.g. hydrogen plasmas |l4ll6lll7lll8l]. helium plasmas [l9|,[2(| and plasmas 
of other elements containing H or He like ions. In most cases, a hot and dense plasma with a large fraction of free 
electrons was considered. In such systems, the main contribution to the shift and width of a given spectral line, 
besides the natural line width and Doppler broadening, comes from the ionic microfield, caused by ions surrounding 
the radiating atom or ion, collisions with free electrons, and interaction with collective excitations in the plasma 
(plasmons). Theoretical results are in very good agreement with experiments carried out, e.g. by Wilhein (2ll. |22||. 

In weakly ionized plasmas, the influence of neutral perturbers becomes dominant. Van der Waals broadening has 
been investigated within second order perturbation theory. Based on the derivation of the van der Waals potential 
between two hydrogen atoms in their ground state as given by London [23| , a red shift of the spectral line follows as 
was shown by Margenau [24| . A review of van der Waals broadening is given by Traving (25j . Line broadening due 
to van der Waals interaction has been investigated in detail, e.g. by Walkup et al. [26j . 

We show how to extend the quantum statistical approach, elaborated for the influence of charged particles on the 
line profile, to include also the interaction with neutral perturbers. In this way, the extension of the van der Waals 
broadening by accounting for many-particle effects such as dynamical screening or strong collisions is possible. A 
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unified description of the influence of the perturbing partially ionized plasma on the line profile of a radiating atom 
can be performed, treating charged particles and neutrals within the same formalism. 

The theory will be compared to optical spectra obtained from hydrogen microdroplets exposed to strong femtosecond 
laser pulses. So far, most of the experimental and theoretical work has concentrated on nanosized clusters [27l. [28l[29| . 
Studies on the light emission from these systems have focused on x-rays and EUV-radiation [3(J HH, [HJ , motivated, 
e.g. by the search for novel light sources. Hydrogen as a target material has attracted considerable attention since 
nuclear fusion has been demonstrated as a result of Coulomb explosion of dense deuterium clusters [33j . Various efforts 
to produce micrometer sized hydrogen targets are reported, cf. e.g. Ref. [34j . Other studies on such large systems 
include, e.g. the measurement of x-ray emission from Kr clusters [35| and from methanol microdroplets [3a. |37I|. An 
expanding microplasma can be produced with interesting parameter values. In particular, strongly coupled plasmas 
are obtained, where density and temperature are time-dependant. Time resolved experimental techniques have been 
demonstrated to measure transient plasma properties like the evolution of the plasma density (38l . |39| . 

In the experiment on /im sized droplets reported here, lines of the hydrogen Balmer series are observed as distinctive 
features in the optical emission spectrum. The present paper will give a first interpretation of the recorded line shapes. 
The analysis of the line profiles shows that van dcr Waals broadening due to neutral perturbers is the predominant 
contribution. Values for the effective plasma parameters are deduced to interprete the observed line profiles, which 
differ for the different lines. However, an equilibrium picture cannot give an agreeable description of the measured 
Balmer spectra. A consistent dynamical picture of the expanding microplasma will be given by use of hydrodynamic 
simulations. This enables us to interpret the observed spectra if considering the time stages where excited atomic 
states can exist. 

The work is organized as follows: In the second section we will briefly outline the method of thermodynamic Green 
functions with emphasis on its application to the calculation of spectral line shapes. The method will be applied to 
the case of a weakly ionized plasma, where the interaction among atoms is governed by the dipole-dipole term. The 
Margenau result appears in the case of ground state perturbers. In the third section, our results will be compared 
to the Balmer spectra obtained from laser produced hydrogen microplasmas. Parameter values for the effective 
temperature and density are inferred. In the fourth section, a hydrodynamical description of the expanding droplet 
after laser excitation is given. 



II. MANY-BODY THEORY OF LINE SHAPES AND VAN DER WAALS BROADENING 



A. Quantum statistical approach to spectral line shapes 

A quantum statistical approach to the optical spectra of dense, strongly coupled systems can be given within linear 
response theory. Emission and absorption of radiation is related to the transverse dielectric function et r (fc, uS). In the 
optical region, the wavelength is large compared with atomic distances so that the long-wavelength limit k — > can 
be considered. For the absorption coefficient we find 

a(u) = ^—lme(0,u) , (1) 
cn(uj) 

n(uj) is the refraction index, c is the velocity of light in vacuum. For k — > the transverse and longitudinal dielectric 
function e(0,o>) become identical. We consider in the following the longitudinal one. 

The dielectric function for a charged particle system can be evaluated in a rigorous way using quantum statistical 
methods. According to e(k,oj) = 1 — (l/e k 2 )H(k,Lu), it is linked to the polarization function H(k,u>), for which a 
systematic perturbation expansion can be given. A Green function approach to the polarization function is outlined in 
App. O Using Feynman diagrams and partial summations, appropriate approximations can be found which account 
for different microscopic processes contributing to the behavior of the charged particle system. In particular, for 
partially ionized plasmas, a cluster decomposition H(k, to) = Hi (A:, to) + Il^ik, u) + . . . can be performed to give the 
contribution of free carriers as well as of bound states in a systematic way [40| . 

Hi(k,uj) collects the single-particle contribution. In lowest order, neglecting collisions with other particles, the 
well-known RPA polarization function 

is obtained. Here, f c (p) = [exp{(i?p — /i c )/fesT} + is the Fermi distribution function of particles of species c 
(including spin) with charge e c and mass m c , E£ = h 2 p 2 /2m c is the single-particle kinetic energy, and /i c the chemical 
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potential. Within this lowest order approximation, single-particle excitations as well as collective plasmon modes are 
described. The limit 77 — > + has to be performed after the summation over momenta. 

The RPA contribution I^ 0) (jfe , u>) is improved if interactions with other particles beyond mean field are included, see 
App.[AJ Medium effects, such as collisions, dynamical screening and exchange, enter the polarization function via the 
single-particle self-energies and the vertex function. The vertex function describes the in-medium coupling of particles 
to the radiation field. It has to be taken in the same approximation as the self-energy as dictated by consistency 
constraints, i.e. Ward identities [4l[. Going beyond the mean-field (Hartree-Fock) approximation, collisions can be 
considered in lowest order Born approximation. Higher orders lead to dynamical screening and t-matrix expressions 
for strong collisions. This way, the continuum of optical spectra is described, in particular inverse bremsstrahlung 
[4^ |. However, line spectra are missing in Hi(k,uj). 

The next term in the cluster decomposition is n 2 (A;,a;). It is given by the convolution of two atomic (two-particle) 
propagators and describes also transitions between bound states, i.e. the line spectrum In the lowest order, any 
interaction with further particles are neglected. Similar to (fc, u>) given above, replacing single-particle propagators 
by atomic propagators one obtains 



l4 0) (fc,c) = 4 



a 1 a 2 ,P 



aiti2 



The factor 4 accounts for spin degeneration, g(E < "°l) — [exp{(E ([ 'i — — ^ e )/fcgT} — 1] 1 is the Bose distribution 
function of electron-ion bound states at E^ p = —Z 2 Ry/n 2 + H 2 P 2 /2M, with a = {n, I, m, m s } denoting the internal 

quantum numbers, HP the center of mass momentum, and M = mj + m c . The unperturbed matrix elements Ma^ a2 (k) 
are given by 
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ip a (p) = (p I a ) and i-'ai^) — (^1 a ) are atomic wave- functions in relative momentum representation and coordinate 
representation, respectively. Eq. (JS]) results for m e jm\ — ► 0. In contrast to a simple chemical picture where a mixture 
of free charge carriers and atoms is considered, scattering states are included, and double counting of diagrams has 
to be avoided. In this approximation, one obtains the unperturbed line-spectrum of isolated atoms including the 
Doppler profile due to thermal motion. 

As before, medium effects are described by taking the self-energy of the constituents and the vertex into account, 
which now are given on the two-particle level, see App.[X] In the simplest approximation, the dynamically screened in- 
teraction can be considered similar to the GW approximation for the single-particle propagator [43[ . This dynamically 
screened interaction contains the polarization function, which once more can be decomposed into the contribution of 
single-particle states, two-particle states and higher cluster states. In this way, we obtain the dynamically screened 
Born approximation for collisions with free particles or composed clusters of the partially ionized plasma. 

In the case of dense and strongly ionized systems, the main contribution to the self-energy comes from the interaction 
of the radiating atom with charged particles, i.e. electrons and ions. The influence of free particles is given by the 
RPA-polarization function, i.e. the first term in the cluster decomposition. This interaction can be inserted into the 
self-energy, but the vertex-function has to be taken in the same approximation. Reducing the screened interaction to 
its one-loop approximation, which amounts to the second Born approximation with respect to the statically screened 
Coulomb potential, the so-called impact approximation is obtained. Systematic improvements by taking higher order 
contributions into account lead to the dynamical screening of the Coulomb potential in the case of weak collisions 
and the t-matrix approach in the case of strong collisions. In particular, the contribution of collisions with electrons 
to the line profile in dense plasmas is treated this way, for a review see Giinter 15 1. 

The following expression for the line-profile £(Aoj) as a function of the frequency displacement Aw = uj—(E^—Ej)/H 
is obtained Hal. 



£(Au;)= £ Mj'(fc) 
H'ff 
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Here, Sj, E/ denote the electronic contribution to the self-energy of the initial or final state, respectively, T v the vertex 
contribution. Improvements of the Born approximation by accounting for dynamical screening and strong collisions, 
have extensively been discussed in the literature [ill ]. 

Due to their large masses, ions are much slower than electrons. In the adiabatic limit, the ions may be considered as 
a static distribution of charged particles during the process of emission. Considering different ionic configurations in 
the plasma, the distribution W p {$) of the so-called microfield (3 — E / En \tsmar\t, -^Holtsmark = (47i ion /15) 2 / 3 Ze/(2e ), 
(n; on is the density of ions with charge Z in the system) is introduced, and the averaging leads to the Stark profile 
as given in Eq. ([5]). Also in the case of ionic contributions, the systematic many-particle treatment leads to further 
improvements such as the dynamic microfield. 

In the next step, one has to take into account the interaction between the radiator and neutral particles (atoms), 
as well. This can be performed along the same lines as for the interaction with charged particles, discussed before. 
The bound state contribution is considered in the cluster decomposition of the polarization function, expanding the 
dynamically screened interaction which enters the self-energy, see App. [A] In particular we will consider the two- 
particle self-energy describing the influence of the medium on the radiator, where the bound state component is 
taken into account in lowest approximation. We consider the second order Born approximation with respect to the 
unscreened Coulomb potential to describe the interaction between neutrals. The detailed calculation is given in the 
next section. 

Again, improvements can be obtained systematically. By calculating the dynamically screened interaction with the 
polarization function up to the second term in the cluster decomposition, the interaction with free and bound particles 
as well as with collective excitations could be included in a consistent way. Strong collisions can be described via the 
four-particle t-matrix. In addition to the self-energy term, also the vertex is modified if bound states are included. 



B. Evaluation of atom-atom interaction in Born approximation 



The calculation of the interaction V aia2 (R) between two atoms at distance R in state \a±) and |aa) is carried out 
in App. [B] We obtain 
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Here, 013, 014 denote quantum numbers of the intermediate states. Only elastic processes are considered where the 
internal quantum numbers for the initial state and the final state are identical. Due to the ion's heavy mass, we 
neglect the contribution of kinetic energy in the intermediate propagator in comparison with the excitation energies. 
Besides R which is the distance between the nuclei of the interacting atoms, f\ and r 2 are vectors of position of the 
two electrons of atom 1 and atom 2, respectively. 

This result ([7]) coincides with results from second order perturbation theory. At large distances Rjas ^S> 1 we can 



perform the dipole approximation using \R — r\ 1 
for R^> r. We obtain the van der Waals potential, 

^4 



l/R 2 R- {R 2 + r 2 - 2R ■ r) /2R + 3(R 2 + r 2 - 2R ■ r) 2 /8i? 2 



E 
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C 2 (^)V'a4(^2) 



(8) 



For hydrogen-like atoms, the sum in Eq. ([8]) has been performed by Eisenschitz and London [23|. Their result for the 
interatomic potential between two hydrogen atoms in ground state (a,\ = a 2 = Is) is given by 



R 6 



(9) 



The constant C§ was calculated as Cq = 12.94 Ry Og, where Ry = 13.6 eV is the Rydberg energy and ae = 0.53 A is 
the Bohr radius. 
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For the interaction between a radiator in an excited state a = n,l and a perturber in ground state, a = Is, the 
following expression for the interaction strength C§{a, Is), can be derived, as shown in App. Id 

o 

C 6 (q,1s)=4 - [5n 2 + 1 - 3Z(Z + 1)] . (10) 

1 + 1 /n z 

Note that there is no exchange contribution in Eq. ([7]), since the distance R is assumed to be large compared with 
the Bohr radius. At shorter distances, the higher order terms as well as exchange terms have to be taken into account 
which will diminuish the van der Waals behavior at small distances and possibly, in the case of electrons in the triplet 
state, result in a repulsion due to the exchange contribution. 



C. Van der Waals profiles 



Within a general approach, the contribution of the dense medium to the spectral line shape is given by the self- 
energy and vertex contribution. The first order of the interaction gives the Hartree-Fock mean field, in second order 
the impact approximation is obtained where collisions are described in Born approximation. 

Within this concept, the van der Waals contribution to the self-energy is 



= E 9(E™p)V£%(q = 0) , (11) 

a 2 ,P 



£vdw = £ „ / d 3 R V£%(R) . (12) 



which in coordinate space reads 



Here, n a2 is the density of pertubing atoms in state ol<x. In the following, we consider the perturbing atoms in their 
ground state «2 = |44j | and only the radiator in an excited state a = nl, n > 1. 

As already discussed in connection with the ionic contribution to the spectral line shapes, because of their large 
masses the atomic motion is slow, and we can apply the microfield concept of strong interaction with a given field 
distribution. The intensity distribution of a given spectral line, due interaction between the radiator in initial and 
final states \i) and | / ), respectively and the perturber in its ground state |0) is obtained by performing an averaging 
procedure over the distribution of perturbing atoms, 



I if (Au) = J dVi Jd 3 r 2 ...J dr N P N (n,r 2 , ...,r N )S ^ Vwfa) ~ Yfofa) - HAuj \ . (13) 

Here, Pn(ti, ...rjsf) is the probability density for having atom 1 in the volume element d 3 ri at r*i, atom 2 in the 
volume element d 3 r 2 at r 2 , etc. The evaluation of the sum (fl"3"|) in the case of statistically independent atoms, i.e. 
P/v(fi, . . .fjsr) — l/V N , is shown in App. IDl 



The result for the dipole limit of the interaction Eq. (J9j) is 24 



77/ w (Aw) = Ai f [-Awf 3/2 e nA *f/ Au} , (14) 

with Aif = (2tt/3) [C^(i, 0) — Ce(f, 0)) 1 / 2 n a t om - This is known as the Margenau profile. Note that Alu is to be taken 
at negative values, i.e. the van der Waals interaction shifts the spectral line to smaller energies (red shift). The 
maximum of the intensity distribution (|14p is located at 

2 7r \ 3 2 

C'e (n a tom a|) Ry/h , (15) 

with Cq = Cq{i, 0) — Cg(/, 0) and n atom the number density of (ground state) hydrogen atoms in the system. The full 
width of half maximum (FWHM) of the line is found as 7 = 2.78 Aw max , see also Refs. [2J, |45| . The Margenau profile 
Eq. ([5]) will be used in order to analyse the hydrogen Balmer spectra obtained from laser produced microplasmas. 
These experiments will be described in the following section. 
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III. EXPERIMENTS WITH LASER-IRRADIATED HYDROGEN DROPLETS 



A. Experimental setup 



Hydrogen droplets are produced by expansion of pre-cooled H2 gas at temperatures between 17 and 22 K through a 
20 /im orifice into a vacuum chamber using a backing pressure of about 8 bar. At these source conditions, H2 droplets 
of about 10 /im in diameter are emitted from the nozzle at a density of rtii q = 4.2 • 10 22 cm~ 3 . Roughly 10 mm behind 
the nozzle the droplets are irradiated by intense IR laser pulses generated by a Ti-Sapphire laser system with a 
wavelength of Ao = 810 nm and a pulse energy of -E pu isc = 2.5 mJ. The beam is focussed to a spot with a beam waist 
of 40 /im. The pulse width r is tuned in the range of r = 50 — 1000 fs, corresponding to laser peak intensities between 
7.0 x 10 15 Wcm -2 and 3.5 x 10 14 Wcm -2 . When the droplets are exposed to the strong laser field, optical emission 
from the hydrogen microplasmas is observed. 
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FIG. 1: Representative emission spectra from laser-excited hydrogen microplasma in two spectral windows for various pulse 
widths r. The spectra are shifted by 1.5 ■ 10 5 counts (a) and 2.5 • 10 5 counts (b) with respect to each other. H Balmer spectral 
lines H a , 656.3 nm, H^, 486.1 nm and H 7 , 434.1 nm are given. The large enhancement at long wavelengths (A > 660 nm) is due 
to Rayleigh scattering of the laser light. 90000 shots are integrated at 1 kHz repetition rate. 



The light is analyzed by a 0.34m grating spectrometer (200 lines/mm) and monitored with a CCD camera, providing 
a spectral dispersion of 0.185nm per pixel. The entrance slit width was optimized for a high signal when integrating 
90000 shots at 1 kHz repetition rate. The resulting spectral resolution of the spectrometer is AA ~ 1 nm. Two regions 
of the visible spectral range are analyzed in more detail, i.e. (a) 430. . .500nm, and (b) 450. . .700nm. Fig. [1] shows 
representative emission spectra from the laser excited hydrogen microplasmas. Elastic Rayleigh scattering of the 
incident radiation produces large enhancement of the measured spectra towards the excitation wavelength Ao . On top 
of the background signal, H-Balmer spectral lines are clearly identified, i.e. H Q at A = 656.3 nm, H^ at A = 486.1 nm, 
and H-y at A = 434.1 nm. 



B. Analysis of the spectra 

1. Effective Temperatures 

In order to quantitatively analyze the Balmer line signals, one has to subtract the background signal. Therefore, 
the neighbouring spectral ranges of each line, i.e. 640. .655 nm and 670. .685 nm for the H Q line, 460. .467 nm and 
472. .480 nm for H/3 as well as 430. .433 nm and 438. .445 nm for H 7 are fitted by an exponential function, which is then 
subtracted from the data set. As an example, Figs. [2] (a) and (b) show the background subtracted spectra for a pulse 
width of r = 1000 fs (dots). The solid curves are the Margenau profiles fitted to the data; they will be discussed in 
more detail in Sec. IIIIB 21 

By comparing integrated intensities from different spectral lines, one can determine the effective temperature in 
the system. Given the ratio of wavelength-integrated intensities I\ / 12 of two spectral lines of the same atomic species 
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FIG. 2: Comparison of experimental line profiles (dots) of H a and H^j from the r = 1000 fs data to Margenau theory Eq. (|14[) . 
taking into account the limited detector resolution A A ~ lnm. For H a , good agreement between theory and experiment was 
obtained setting n^ m = 4.5 • 10 20 cm" 3 . For H,g, the best fit was obtained with n^ m = 1.8 ■ 10 20 cm -3 . 



and the same ionization stage at central wavelengths Ai and A2, and assuming again local thermal equilibrium, the 
effective temperature is [45| 



J 12 — 



E± — E% 



fee In 



A% 31 A 2 L2 



A 2 gi Ai h 



(16) 



Ai, A2 are the angular momentum averaged Einstein coefficients for the transition between levels corresponding to 
the spectral line 1 and 2, respectively and g\ and (72 are the degeneracy factors of the excited level (principal quantum 
number nf x ) of each transition, i.e. gi — 2(n° x ) 2 . E\ and E2 are the energies of the excited level for each transition. 
Tab. Q] lists effective temperatures for each pulse width. As can be seen, is dependent on r. One obtains different 
results from the two ratios I(H a )/I(lIp) and 7(H ( 3)//(H 7 ). For small r, rises with increasing r. At t ~ 400 fs, 



T°j| saturates near 1 eV. 



On the other hand, TJ? gives values around 0.5 eV. 



TABLE I: Effective line temperatures according to Eq. (|16p versus pulse length. The first set of figures was obtained from the 
short wavelength spectrum, shown in Fig. QJa), the second set results from the data shown in Fig. [TJb) . 



r[fs] 


60 103 233 


r[fs] 


50 109 170 264 392 695 1000 


Tf, [eV] 


0.49 0.56 0.41 


Tf B [eV] 


0.57 0.62 0.81 0.97 0.99 1.00 1.04 



These effective temperatures have to be analyzed within a microscopic description of the expanding plasma source. 
Self-absorption should be considered so that the intensities occurring in Eq. (fTS)) are modified when propagating 
through the microplasma. Furthermore, different values from the intensities of H a , Hg and H 7 indicates that the 
system is in a generic non-equilibrium state, which cannot be described by a single temperature. Despite the measured 
spectra are time integrated, we have to assume a dynamic evolution of the plasma parameters temperature and density 
over the course of the expansion of the excited hydrogen droplets. Within a hydrodynamical description given below 
in Sec. IIV AI local thermal equilibrium with plasma parameters depending on space and time are introduced. 

During the expansion process, radiation can only be produced if the upper level of the transition in question is 
existent in the system so that it can be occupied. From simple geometric and energetic arguments it is clear, that 
n = 3,4,5 levels are not present in hydrogen at liquid densities. This point will be discussed in more detail in 
Sec. IrTBl 

The continuum background stems from free-bound and free- free (bremsstrahlung) transitions. The spectral behavior 
can also be used to estimate effective temperatures which should be analyzed within a dynamical expansion model. 
Since more experimental work is necessary to separate the continuum background radiation, we will consider only the 
line spectra in this work. 
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2. Spectral line shapes 

To analyze the line profiles, one first has to consider the effect of the broadening due to the finite resolution of the 
experimental setup. The true line profile is obtained after deconvolution with the detector function. In a separate 
experiment the spectral resolution was measured as AA ~ 1 nm . This is consistent with the steep increase of intensity 
on the blue wing of the unperturbed spectral line. 




656 658 660 662 664 666 



wavelength A, [nm] 

FIG. 3: Normalized spectra near the H a -line for r = 1000 fs pulse length. The raw data (cf. inset) have been summed over 
different vertical ranges. The solid line marks the central region of the laser focus (framed part in the inset) , the dashed curve 
corresponds to the data from the peripheric region. 

On the long wavelength wing, all spectral lines show a significant broadenening. In Fig. [3J the spectral line shape 
is analyzed as a function of the distance of the origin of radiation from the plasma's center. This is achieved by 
integration of the CCD-Data over different ranges in the vertical direction, cf. the inset in Fig. [3J The corresponding 
spectra for the central region of the laser focus (framed part in the inset) and for the peripheric region, about 100 /im 
beside the focus are shown. After background subtraction, the spectra are normalized to their maximum value. In 
this way it can be seen, that the width of the spectral line depends significantly on the spatial origin of the radiation 
and decreases with increasing distance from the center. This behaviour indicates that the asymmetric broadening of 
the spectral lines is related to the properties of the microplasma and should be described as a density effect. 

Pressure broadening of spectral lines is caused by charged as well as by neutral perturbers. As discussed in Sec. HT1 
free electrons are treated in impact approximation leading to a symmetric Lorentzian line profile. Also the ionic 
microfield contributes to a Voigt profile with symmetric broadening on both the red as well as the blue wing (linear 
Stark effect). Both effects are smaller than the spectral resolution of 1 nm, as can be seen in the spectra, e.g. Fig. [3J 
Thus, we conclude that the influence of free charged particles is not clearly identified so that the free electron density 
is below n f l cc = 10 17 cm" 3 which at T = 1 eV gives for H Q the FWHM of 1 nm. 

The asymmetric red shift can be described by the interaction with neutral perturbers. We will now compare the 
experimental data to the line-shape due to interaction with neutral perturbers, as outlined in Sec. Ill Bl In Fig. [21 the 
measured spectra in the vicinity of the H Q line are compared to the emission profile as given by Margenau, Eq. (|14p . 
convoluted with the detectors resolution function In the case of the H a -line, the best fit was obtained using 

n itom = 4-5 • 10 20 cm~ 3 for the density of neutrals. In the case of Hp, the density nf2, m — 1-8 ■ 10 20 cm" 3 gives the 
minimum \ 2 - 

As for the effective temperatures, also the different values for the effective density contradict an equilibrium picture. 
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A dynamical description of the microplasma's evolution is needed. This task will be accomplished in the next section 
by means of hydrodynamical simulations. 

IV. DYNAMICAL PLASMA EXPANSION MODEL 

In the previous section the effective temperature and density of the microplasma were determined by analysis of the 
experimental spectra. For both quantities, different values have been obtained, depending on the external parameters, 
e.g. the pulse width of the laser, and on the wavelengths considered in the analysis. This strongly indicates that the 
plasma parameters are time-dependant quantities, which have to be modeled, e.g. by hydrodynamic simulations to 
improve the description of the microplasma and the interpretation of the experimental data. 

A. Hydrodynamic expansion 

Hydrodynamic simulations are a versatile tool to infer the dynamics of a strongly coupled many-particle system. In 
the case of plasmas, hydrocodes have been successfully applied to study these systems under the influence of strong 
external fields [47j as well as the relaxation of a plasma in an excited state into equilibrium. Here, we will study 
the hydrodynamic expansion of the excited, i.e. heated H microdroplet after times which are long compared to the 
pulse length of the laser. Simulations are performed using the hydrocode MULTI2002 [H, |49j|. As initial conditions, 
we assume a homogenous density profile n(r, t = 0) = ny lq = 4.5 • 10 22 cm~ 3 and also homogeneous temperature 
distribution. Assuming 50 % absorption of the laser energy given above by the droplet [501 ] and accounting for the 
ratio of the droplet's cross-section to the focal spotsize, i.e. (10 /im/40 /im) 2 ~ 0.06 as well as for the energy needed 
to break up the molecular bounds, i.e. 2.2 eV per atom, we obtain T(r, t = 0) = 14 eV as the initial temperature. 




radius [urn] radius [jun] 

FIG. 4: Hydrodynamic simulation of the density (left) and temperature (right) profiles of an expanding hydrogen microplasma 
starting from nu q = 4.5 • 10 22 cm -3 and T = 14 eV. The broad solid and dashed curves give the conditions after 0.7 ns and 0.9 
ns. At this time, the central density has reached the density which has been obtained from the fit of the Margenau profile to 
the H a and the spectral line, respectively for r = 1000 fs pulselength. The corresponding temperatures at these times are 
close to the effective temperatures T°j| ~ 1 eV and TI? ~ 0.5 eV obtained from the spectra, see Tab. HI 

Fig. U shows the number density and temperature profiles. The uppermost curve corresponds to t = ns, while the 
last curve is taken at t = Ins. Both parameters n(r,t) and T(r,t) show a fast decline in the center of the droplet. 
The density decreases by a factor of 10 every 300 ps. Secondly, all profiles are nearly constant over several tens of 
micrometers starting from the center and decrease sharply when approaching the rim of the droplet. Therefore, the 
central values n(r — 0,t) and T(r = 0,t) can be taken as representative values for the whole droplet. The temporal 
evolution of the central number density and temperature is shown in Fig. [51 
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timet[ns] time t [ns] 



FIG. 5: Temporal evolution of central density n(r = 0, t) (left) and temperature T(r — 0,t) (right), given as bold solid lines. 
Narrow solid line: central parameters after 0.7 ns, i.e. at the time, when the density obtained in the fit of the Margenau profile 
to the H a line, = 4.5 • 10 20 cm -3 is reached. Dashed line: the same for Hp. 




time t [ns] 

FIG. 6: Ionization degree a e as a function of time obtained from Saha equation (bold solid curve). The narrow horizontal line 
gives the ionization after 0.7 ns, i.e. at the time, when the central density has decreased to = 4.5 ■ 10 20 cm -3 , which has 

been obtained as best-fit parameter with respect to the Margenau profile of the H a spectral line. 

The hydrodynamical description gives the time evolution of the excited droplet. The assumption of local thermal 
equilibrium may be justified on the time scale of ps for the thermalization of kinetic energies. If assuming also local 
ionization equilibrium, the composition during the time evolution can be calculated. The temporal behaviour of the 
ionization in the droplet's center is plotted in Fig. [5] For the initial conditions T = 14 eV and nu q = 4.2 • 10 22 cm -3 , 
we have an ionization degree a e = ri{ Tee /n = 0.68. After t = 0.7 ns, the density and temperature have decreased so 
far, that a e drops below 10~ 4 . This corresponds to a concentration of free electrons of roughly 10 17 cm~ 3 . At this 
density, and T ~ 1 eV, the Stark effect leads to a broadening of the H a line of 1 nm, as was discussed in Sec. IIIIB 21 
Thus, at times larger than 0.7 ns, Stark broadening does not give a notable contribution to the width of the spectral 
lines. 

This analysis has shown, that the emission of Balmer line radiation occurs at comparatively late times of the 
expansion, when the density has decreased to a certain level. This is due to the fact, that in a dense system, excited 
levels are not well defined due to interaction with neighbouring atoms. Both wavefunctions and atomic potentials are 
disturbed. In the following section we will analyse the question, which density has to be established in the system, so 
that excited atomic levels are defined and the corresponding radiative transition may occur. 
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B. Occupation of excited levels 



In a dense system, the potential energy of an electron in a given atom is modified by the medium [51(. Thus, energy 
eigenvalues and wave functions are changed due to screening by free charge carriers. In particular, it is well-known 
that bound states merge with the continuum of scattering states at high densities so that the electrons are no longer 
bound to a special ion, but move relatively free in the plasma as denoted by a transition from dielectric to metallic 
behavior [H^] • Considering excited states, the dissolution occurs already at lower densities so that these excited states 
cannot any longer be occupied to produce a line spectrum. 

We focus here to the influence of bound states in the medium to analyze at which densities the excited states are 
dissolved into the continuum so that no line spectra are formed. In a first approximation, we calculate the modification 
of the potential in the Schrodinger equation for the hydrogen atom due to the perturbing atoms in the surrounding 
medium. Only ground state atoms are considered, and the correlations between the electrons in the perturbing 
ground state with the radiating electron in the excited state are neglected. Solving the Poisson equation, the effective 
potential is obtained as superposition of the Coulomb potential and the potential of the neighboured atoms, which 
leads to a lowering of the Coulomb potential by the amount 




AVi(f) = f!__ e -2|^l/«B | | . r l^i j , | 7) 

47reo|f— R4 

Here, R4 is the position of the perturbing atom. The lowering of the Coulomb potential has two effects: On the one 
hand, the energy level is shifted to lower energies, on the other hand, at a certain critical distance, the threshold of 
potential energy between neighboured atoms crosses the shifted electrons energy level, and the bound state is dissolved 
into the continuum. This is the case, if the distance between neighbours comes into the range of the extension of 
the electrons wave-function. We assume a closely packed configuration where the radiating atom is surrounded by 12 
perturbing next neighbours. The shift of the energy levels with principal quantum numbers n = 3 and 4 has been 
evaluated in first order of perturbation theory and can be given in an analytical form. Results for the shift as a 
function of the interatomic distance are given in Fig.[7]Ja). 




FIG. 7: (a) Energy shift of 3s and 4s levels as a function of the distance R to the next neighbouring atom, 
(b) Effective potential due to the overlap between the radiator's Coulomb potential and the screened potential of the next 
neighbouring atom in its ground state. Results are shown for two different values of the interatomic distance R = 20(Jb 
(bold curve) and R = 30 Ob (dashed curve). The perturbed 3s (R = 20 Ob) and 4s (R = 30 Ob) energies are given as narrow lines. 



In Fig [7jb) we give the potential energy V(r) in the direction to a next neighbour for two different distances, 
together with the corresponding shift of the energy levels with n = 3 and n = 4. As can be seen, the bound state 
disappears at a critical distance R^ where the threshold energy becomes lower than the binding energy, and the 
electron is no longer bound to the central ion, but escapes the effective potential and moves relatively freely within 
the cluster. 
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From the critical distance Rc , one can then infer the "critical" density ni n> that has to be established in the 
expanding system, before the upper level n of the considered transition can exist, i.e. before the corresponding spectral 
line can occur. In this way, we obtained for the second excited level, relevant for the H Q transition, a "critical" density 
of ni a) ~ 1.24 -10 21 cm 3 . This value for n c is by a factor of two larger than the value n^l m = 4.5 • 10 20 cm 3 obtained 
from the fit of the H a line to the line shape due to van der Waals interaction, cf. Fig. [2ja). 

It has to emphasized in this context, that these are exploratory calculations to estimate the region of density where 
well-defined excited energy levels of the radiating atom can exist. A more detailed calculation should include electronic 
correlations between radiator and perturber which have been neglected in the calculation of the perturbing potential. 
The account for correlation and exchange effects will modify the potential and the critical density. In addition, we 
considered only a perturbing atom at mean distance, neglecting any fluctuations in configurations. On the other hand, 
the experimental spectra are time integrated measurements and hence the density has to be interpreted as a mean 
value, averaged over the whole exposure time. At the moment the radiation starts, the density might in fact be larger 
than the value inferred from the measured spectra. For Hp, the critical density is obtained as ni^ ~ 2.66 • 10 20 cm~ 3 , 

which also exceeds the fit parameter for the line (n[2 m = 1.8 ■ 10 20 cm~ 3 ) for the same reasons as given above for 
the case of the K a line (cf. Fig. 0(b)). 

In conclusion, the experimental Balmer spectra are not characteristic for the first stage of the laser excitation of 
the cluster, because the density of hydrogen in the condensed state is too high to form well-defined excited atomic 
levels with the corresponding principal quantum numbers. Such levels appear only during the process of expansion 
and may be used as a signal to infer the state of the microplasma at the corresponding time stage. 

C. Line emission scenario 

As discussed in the previous Section, the observed Balmer spectra cannot be interpreted within an equilibrium 
picture of the laser produced microplasma. We could not infer consistent values of plasma parameters for temperature 
and density. A consistent description is only possible if the time evolution of the expanding microplasma is considered 
as a non-equilibrium process. 

We follow the dynamical expansion of the laser excited hydrogen droplet as given by the hydrodynamical calculation. 
First, the time evolution of the density is considered, with account to the effective densities derived from the measured 
line profiles. Thereby, one can determine the time at which the H Q line appears, i.e. where the inferred density 
n itom — 4-5 • 10 20 cm~ 3 is reached. Starting with the excitation temperature T — 14 eV after the short-pulse laser 
excitation, this density is established after 0.7 ns. The corresponding temperature obtained from the hydrocode 
MULTI is T = 0.65 eV and is in reasonable agreement with the estimate for T®5 given in Sec. IIIIB 21 Thus, the H a 
line profile reflects the state of the expanding microplasma after 0.7 ns. In the case of the H^ we observe the following: 

From the van der Waals fit to the data, cf. Fig. [Hb), = 1.8 ■ 10 20 cm -3 was obtained, which is reached in the 

hydrodynamic simulation after t ~ 0.9 ns. The temperature is T ~ 0.4 eV at that moment, which is close to the value 
0.5 eV obtained for TJ^, cf. Tab.D thus gives us information about the state of the droplet at t ~ 0.9 ns. 

Second, as for the empirically determined plasma parameters, the critical densities for H Q and as discussed in 
Sec. IIVB1 have to be compared to the results of the hydrodynamic simulation. The central density of the droplet 
decreases to the critical density for H Q to appear, ric = 1.24 ■ 10 21 cm' 3 after 0.45 ns. After the same time, the 
temperature at the center of the droplet has decreased to 1.3 eV, which is in good agreement with the value 1 eV 
inferred from T^jg. The critical density for H3, ni^ = 2.66- 10 20 cm -3 is reached after 0.8 ns. The central temperature 
is 0.5 eV at that moment, which again coincides well with the value 0.5 eV as obtained from cf. Tab.|TJ 

This analysis shows, that the droplet in fact undergoes a complex dynamical evolution. Using hydrocodes to 
model the droplets history, we can understand the experimental observations. The initial temperature of 14 eV is 
consistent with the observed temperatures and densities established at later times of the evolution. The line radiation 
obviously stems from relatively late times. What happens before line radiation appears? At the initial conditions 
of riatom = n\iq = 4.2 • 10 22 cm~ 3 and T — 14 eV, the ionization of the droplet was calculated as a e — 68%, solving 
Saha's equation, cf. Fig. [6] The plasma emission is thus dominated by continuum radiation (free-free and free-bound 
transitions). However, since the detector integrates over the whole evolution of the plasma, only a nearly constant 
background remains in the spectrum from this early stage of the evolution. As the ionization decreases, bound 
states begin to form and line radiation occurs. However, for a given spectral line, both upper and lower level of the 
corresponding transition have to be well defined. To this end, the density has to fall below a certain value, in order 
to allow for the excited energy level to appear underneath the effective potential in the dense medium. 

Considering the time scales for the expansion of the cluster of the order 10 -9 s according to the hydrodynamical 
calculations, the local thermodynamic equilibrium can be assumed, and the life time of the excited states of the 
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order 10 11 s is short compared with the evolution of the microplasma. In a more rigorous approach, the variation of 
temperature and density in space and time should be accounted for to synthesize the spectra. 



V. CONCLUSIONS 

We found that the visible spectra can be used to get signatures for a time-resolved picture of an expanding mi- 
croplasma. The line profiles allow for the determination of the microplasma's temperature and density during its 
hydrodynamic expansion. We considered an energy deposition in the liquid hydrogen droplet which is sufficiently 
weak so that in the expanded phase most of the electrons are found in atomic bound states. Thus, the spectral line 
profiles are determined by van der Waals broadening. We presented a general quantum statistical approach to line 
profiles which allows for the unified description of charged and neutral perturbers, including density effects such as 
dynamical screening or strong collisions. 

Although the spectra are time integrated measurements, we can use the Balmer lines as signatures for the expanding 
microplasma at definite values of density. A consistent scenario for the emission of Balmer lines has been given. To 
follow the time evolution of an excited droplet more directly, time resolved spectra have to be analyzed. Possibly, this 
can be achieved with pump-probe experiments and/or stimulated emission measurements. 
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APPENDIX A: GREEN FUNCTION APPROACH TO THE POLARIZATION FUNCTION 

The general expression for the single state contribution IIi(fc, uj) to the polarization function has the form 



G 1 




n!(fc,w M )= cf >• (Ai) 



The full single-particle propagator contains the self-energy Si(p, z), 



+ — ' W*— . (A2) 



Gi G< 0) G< 0) E x G, 

The vertex V describes the coupling to the electromagnetic field and can also be expressed in terms of an effective 
interaction kernel. Both quantities, self-energy X and vertex T, have to be approximated in a consistent way. For the 
self-energy one has the GWF" approximation, given by the diagram 




{p,z v )= e X , (A3) 
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which is a self-consistent equation for the full propagator, the full vertex T, as well as for the screened interaction 
potential W, 




V 




n >v w 



V{k) 



i-y(fc)n(fc,^) 



(A4) 



while the vertex function Y is the solution of the Bethe-Salpeter equation, given in diagrammatic form by 





(A5) 



Solving Eq. (|A3[) and Eq. (|A5|) simultaneously is a formidable task. Considerable simplification of the problem is 
obtained by replacing the full vertex V in Eq. (|A3[) by the bare vertex T^°\ This is the so-called GW approximation, 





(A6) 



r (o) G 



The GW approximation is well studied in condensed matter physics. It also describes bremsstrahlung an can be 
further improved accounting for plasma effects, see Ref. (43|. 

The two-particle contribution n2(fc,w) to the polarization function is described in analogy to Eq. (|A1|) . 



Go 



n 2 (fc, w M ) 



where we have the full two-particle propagator 




G 2 
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(A7) 



G 2 {aP,n x ) = 
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with the two-particle self-energy T, 2 (aP, 




(A8) 



(A9) 
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and the vertex follows from an equation in analogy to Eq. (IA5|1 




(A10) 



For the polarization function which defines the screened interaction used in the two-particle self-energy Eq. (|A9[) . we 
perform the cluster decomposition as outlined in Sec. [TXJ In the impact-approximation, we replace the full screened 
interaction by the first non- ideal term in the iteration of Eq. (|A4[) . Diagrammatically, the two-particle self-energy in 
impact-approximation reads 



s imp K,o A ) = 




(All) 



Note that one has to take into account that double counting has to be avoided so that scattering states in H2(k,oj) 
should not interfer with contributions from IIi(fc, ui). Important are the account of bound states in Il2(k,uj). 

The higher order corrections to the full vertex M(q) can be shown to effectively reduce the two-particle self-energy 
(lH by the amount iY" . After analytic continuation u)^ — > ui + i0 + to real frequencies, we obtain 



Q1Q2 ,P 



uj- E 



(0) 



£ (0) _-Rc 



a 2 P-k aiP 



(A12) 



Note, that the dependance on ui is neglected in the self-energy of the lower energy level [T5|. Neglecting the self-energy 
and the effective vertex in Eq. (|A12j) . we arrive at Eq. ([3]). 



APPENDIX B: BOUND STATE CONTRIBUTION TO THE TWO-PARTICLE SELF-ENERGY 



The interaction between two atoms in state |ai) and \ai), moving with momentum Pi and P 2 respectively, and 
carrying Matsubara frequencies Q\ and f2 K , upon exchange of momentum k is given by the following diagram: 



AfW Gf ] M<°> G 2 0) 



V i V 

I 

M (o) G (0) M ( 0) G (0) 



£ M «L (fc - fc')M(° 2 L 4 {-k')Mi% 2 (k + k')V n , V n _ p 



A + ft K - E 



(o) 



(0) 



(Bl) 



(B2) 



with the unscreened Coulomb propagator Vg = — l/eofc 2 Slo- Due to the ion's heavy mass, we neglect the transfer 
momentum k in the intermediate propagator as well as the kinetic energy h 2 P 2 /2M. Replacing the Matsubara 
frequencies Q\ and f2 K by their on-shell values, we can perform the summations over momenta in Eq. ^ and arrive 
at the familiar expression Eq. (7|). 



17 



APPENDIX C: CALCULATION OF THE INTERACTION STRENGTH 

The expression Eq. {§]) has to be evaluated. We rewrite it as 



yvdw /m _ C 6 (ai,a 2 ) 

v a 1 ,a 2 \ JX ) jj>6 

.4 



E' 



|(aia 2 |W| a i a ii 



,2 (CI) 



(0) 



(47re ) 2 i? 6 ^ K (0) + E {0) - K (0) - E 

The primed sum indicates summation only over states, which assure a finite denominator. 

Henceforth, we assume one atom in an excited state ai = n\l\, ri\ = 2, 3, 4, ... , and the second atom in its ground 
state, a 2 = 0. The denominator is dominated by the term Ea{ + Ea^ = — Ry(l + l/nf). Eq. (|Cip becomes 

/ 2 \ 2 1 



^° )= fej (l + l/^)R y E'K^|Wk3a 4 )| 2 (C2) 
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^ (aid |W| a 3 a 4 )(a 3 a4 |W| aiO) - |(ai0 |W| ai0)| 2 (C3) 
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The diagonal matrix element of W vanishes for atoms with no permanent electric dipole moment. Chosing the 
nucleus-nucleus axis parallel to the ir-axis of the coordinate system, the operator W reads 

W = yiy 2 + ziz 2 - 2xix 2 . (C5) 

Since the expectation values of any coordinate or z vanishes, e.g. (a%0 |ri| a\0) = 0, i = 1, 2, 3, the mixed terms 
in W 2 vanish. We obtain 
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V47reo7 (l + l/7i 2 )Ry 
(0 |y 2 | 0)<ai lyilax) + (0 |z 2 | 0)<ai |z 2 | Ql ) + 4(0 |x 2 | 0)(ai |x 2 | a x ) , (C8) 



and with 



(rf) = J(r 2 ) (C9) 



Ce(ai ' 0) = (i^) (iTT7^% l <0|r2|0)<Ql|r ^ |ai) ■ (C10) 
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Using 



2 2 

(nl\i 2 \nl) = [5/i 2 + 1 - 3Z(Z + 1)] , (Cll) 



we finally obtain Eq. (fT0|) . 

Note, that for the case, that the first atom is also in its ground state, i.e. a = 0, Ce(0,0) = 12 is obtained. The 
difference to the London-Eisenschitz result 12.94 is due to the approximative treatment of the denominator in Eq. (jClj) 
in the calculation presented above. For excited states, which are considered in this work, this contribution becomes 
negligible. 



APPENDIX D: EVALUATION OF THE VAN DER WAALS PROFILE 



We start with the general form of expression (O, i.e. the intensity distribution of a given spectral line, due to 
interaction between the radiator, which is in state | i ) before the transition and in state | / ) after the transition and 
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the perturber i, i = 1 . . . N in state |rij) . It is obtained by performing an averaging procedure over the distribution 
of the N perturbing atoms, 



Iif(Aw) = J dVi J d 3 r 2 ...J dr N P N ( 



N 

r N )5 ( ^2V iaj (rj) - V faj (rj) - Acj 

3=1 



(Dl) 



Here, P/v(fi, . . . rj\r) is the probability density for having atom 1 in the volume element d 3 fi at fx, atom 2 in the 
volume element d 3 r2 at r*2, and so forth. In can be expanded in a cluster decomposition as 



fW(fi,f2, •••>fv 



\ j<k j<k<l 



ffc.n) + 



(D2) 



where the lowest non-ideal term is the pair distribution function, giving the probability to find a second particle at 
fk if there is one particle at f / . 

The 7V 3 -fold integral in Eq. (|D1|) can be evaluated using the Fourier representation of the delta function, 



1 f°° 

8($ - Auj) = — I dpi 

27r J-oo 



ip($ — Auj) 



(D3) 



with <!> = Y] _j (rj) — Vf aj [rj)]. In the case of statistically independant atoms, which are all in the same quantum 



state a.j = a p i.e. Pjv(fi,f2, ■ ■ ■ jf/v) — V , one immediatly finds 
Using the van der Waals limit of the atom-atom potential, Eq. ([9|), Eq. (|D4j) turns to 



dpe 



-iAujp 



N 



dpe 



-iAup 



drr 2 e- lpC6/r6 



N 



The integration over r can be rewritten in the form 
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(D4) 

(D5) 

(D6) 
(D7) 



with Cq = Cq(/,0) — C§{i,0). In the thermodynamic limit, i.e. V — > oo,N — > oo, while keeping constant the density 
of neutrals iV/V = ?i a tom = const, one obtains 
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For V(p), 



is obtained, and finally the Fourier integral Eq. (|D5|) can be evaluated to give 

27T^atomV M 

3(-Aw)3/ 2 



27 TO atom ^fCl ( 27r „ atom / 3 ) 2 / Aa 



(D9) 
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which corresponds to Eq. (fl4|) . 



